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A number of basic scientific questions relating to ion conduction in homogeneously disordered 
solids are discussed. The questions deal with how to define the mobile ion density, what can be 
learned from electrode effects, what is the ion transport mechanism, the role of dimensionality, 
and what are the origins of the mixed-alkali effect, of time-temperature superposition, and of the 
nearly-constant loss. Answers are suggested to some of these questions, but the main purpose of the 
paper is to draw attention to the fact that this field of research still presents several fundamental 
challenges. 
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I. INTRODUCTION 

Ion conduction in glasses, polymers, nanocomposites, highly defective crystals, and other disordered solids plays an 
increasingly important role in technology. Considerable progress has been made recently, for instance with solid-oxide 
fuel cells, electrochemical sensors, thin- film solid electrol ytes in batter ies and supercapacitors, electrochromic windows, 
oxygen-separation membranes, functional polymers, etc! 1 * 2 * 3 * 4 * 5 * 6 * 7 * 8 * 9 * The applied perspective is an important catalyst 
for work in this field. In this paper, however, the focus is on basic scientific questions. This is relevant because ion 
transport in disordered materials remains poorly understood. There is no simple, broadly accepted model; it is 
not even clear whether any generally applicable, simple model exists. Given the intense current interest in the 
field - with hundreds of papers published each year - it is striking that there is no general consensus on several 
fundamental questions. 10 This is in marked contrast to other instances of electrical conduction in condensed matter 
where a much better understanding has been achieved, e.g., for electronic conduction in metals, semiconductors, and 
superconductors, as well as for ion conduction by defects in crystals. 

This paper su mmarizes and discusses bas ic scientific questions relating to ion conduction in (mainly) homogeneously 
disordered solids^ ^ 2 * 13 * 14 * 15 * 16 * 17 * 18 * 19 * 20 * 21 *. The main motivation is not to suggest or provide answers, but to inspire 
to further research into the fundamentals of ion conduction in disordered solids. A question that is not addressed 
below, which has been a point of contro versy pa rticularly during the last decade, is how to best represent ac data, via 
the conductivity or the electric By now this has been thoroughly discussed in the literature, and 

we refer the interested reader to the discussion in Refs. 26 27 28] that present and summarize the differing viewpoints. 



II. HOW TO DEFINE MOBILE ION DENSITY? 

Ion motion in disordered solids is fundamentally different from electronic conduction in crystalline solids. Ions 
are much heavier than electrons so their motion is far less governed by quantum mechanical descriptions. Below 
typical vibrational frequencies (< 100 GHz) ion motion can be described by activated hopping between (usually) 
charge compensating sites. Moving ions carry charge, of course, and thus produce an electrical response which can 
be detected by a variety of experimental techniques. Unlike crystals, the potential energy landscape experienced by 
an ion in a glass or otherwise disordered solid is irregular and contains a distribution of effective depths and barrier 
heights, as sketched in Fig. 1. The effective energies result from differing binding energies at residence sites and 
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FIG. 1: Schematic figures illustrating ion jumps in a disordered landscape, here in one dimension. The arrows indicate attempted 
jumps. Most of these are unsuccessful and the ion ends back in the minimum it tried to leave: If the barrier is denoted by 
AE, if T is the temperature and ks is Boltzmann's constant, according to rate theory the probability of a successful jump is 
roughly exp(— AE /ksT). This implies that on short time scales only the smallest barriers are surmounted. As time passes, 
higher and higher barriers are surmounted, and eventually the highest barriers are overcome. In more than one dimension the 
highest barrier to be overcome for dc conduction is determined by percolation theory; there are even higher barriers, but these 
are irrelevant because the ions go around them. 

differing saddle point energies between residence sites, and they are influenced by interactions between the ions. With 
increasing time scale, the ions can explore larger parts of space by overcoming higher energy barriers. 

Following standard arguments, suppose ions with charge q are subjected to an electric field E. The field exerts the 
force qE on each ion, resulting in an average drift velocity v in the field direction. The ion mobility \i is defined by 
fi = v/E. If the number of mobile ions per unit volume is n mo b, the current density J is given by J = gn mo b^- Thus 
we obtain the following expression for the dc conductivity defined by <j^ c = J/E: 

&dc qrimob/J, . (1) 

This equation expresses the simple fact that the conductivity is proportional to the ion charge, to the number of 
mobile ions, and to how easily an ion is moved through the solid. As such, Eq. ([!]) is an excellent starting point for 
discussing how the conductivity depends on factors like temperature and chemical composition. Or is it? We shall 
now argue that the above conventional splitting of the conductivity into a product of mobility and mobile ion density 
involves non- trivial assumptions. 

Except at very high temperatures ion motion in solids proceeds via jumps between different ion sites. Most of 
the time an ion vibrates in a potential-energy minimum defined by the surrounding matrix. This motion does not 
contribute to the frequency-dependent conductivity except at frequencies above the GHz range; only ion jumps 
between different minima matter. The mobility reflects the long-time average ion displacement after many jumps. 
The fact that ions spend most of their time vibrating in potential energy minima, however, makes the definition of 
mobile ion density less obvious: How to define the number of mobile ions when all ions are immobile most of the 
time? 

Intuitively, Eq. still makes sense. Imagine a situation where some ions are very tightly bound ("trapped") while 
others are quite mobile. In this situation one would obviously say that the density of mobile ions is lower than the 
total ion density. The problem, however, is that the tightly bound ions sooner or later become mobile, and the mobile 
ions sooner or later become trapped: By ergodicity, in the long run all ions of a given type must contribute equally to 
the conductivity. Thus on long time scales the "mobile" ion density must be the total ion concentration. This "long 
run" may be years or more, and ions trapped for so long are for all practical purposes immobile. Nevertheless, unless 
there are infinite barriers in the solid, which is unphysical, in the very long run all ions are equivalent. 

The question how many ions contribute to the conductivity makes good sense, however, if one specifies a time scale. 
Thus for a given time r it is perfectly well-defined to ask: On average, how many ions move beyond pure vibration 
within a time window of length r? If the average concentration of ions moving over time r is denoted by n mo b(T) and 
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FIG. 2: (a) Schematic figure showing the real part of the ac conductivity as a function of frequency at three different temper- 
atures. As temperature is lowered, the dc conductivity decreases rapidly. At the same time the frequency marking onset of 
ac conductivity also increases (actually in proportion to the dc conductivity, compare to the BNN relation Eq. (13) discussed 
below), (b) The real part of the conductivity at three different temperatures for a Lithium-Phosphate glass. 



n is the total ion concentration, ergodicity is expressed by the relation 



^mobO" ^ °°) 



(2) 



An obvious question is how to determine n mo b(r) experimentally. A popular method of determining the "mobile 
ion density" - without explicit reference to time scale - is by application of the Almond- West (AW) formalis m 129 * 30 * 31 * 
that takes advantage of the frequency dependence of the conductivity. We proceed to discuss this approach. First 
note that in ion conductors with structural disorder, the short-time ion dynamics is characterized by back-and-forth 
motion over limited ranges, "sub diffusive" dynamics, whereas the long-tim e dynamics is characterized by random 
walks resulting in long-range ion transport, "diffusive" dynamics (Fig. l)) 32 l 33 l 34 l 35 l 36 l The back-and-forth motion 
leads to dispersive conductivity at high frequencies, while the long-range transport leads to the low-frequency plateau 
marking the dc conductivity (Fig. 2). There is experimental evidence that in materia ls wit h high ion concentration, 
at any given time only part of the ions are actively involved in back-and-forth mot ion .Ell 3 ^. 

A widely applied description of conductivity spectra in the low- frequency regime (i.e., below 100 MHz) is a Jonscher 
type power law, 



(3) 

where we have written the equation in a form such that the crossover frequency marking the onset of ac conduction, 
i/*, is given by a f (i/*) — 2<Jd c - Equation (3) is sometimes referred to as the Almond and West (AW) formula, although 
Almond and West did not consider Jonscher's "universal dielectric response" of disordered systems, but introduced 
their formula to describe defective crystals with an activated number of charge carriers. Nevertheless, when applying 
Eq. (3) to strongly disordered systems, as, e.g. ionic glasses, many authors in the literature follow the physical 
interpretation suggested by Almond and West and identify the crossover frequency with a "hopping rate". Thus 
combining this ansatz with the Nernst-Einstein relation gives 



6k B T a dc 

%w = — — (4) 
q z a z v* 

as an equation to determine the number density of "mobile ions", uaw (after typically assuming jump lengths 
a = 2 - 3 A). 

However, if one accepts that Eq. (|3| provides a good fit to spectra in the low-frequency regime - it generally fails 
at frequencies above 100 MHz - the estimate of an effective number density of "mobile ions" based on Eq. Q is 
questionable. Application of the fluctuation-dissipation theorem implies the following expression, where t* = and 
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H and 7 are numbers that are roughly of order unity (H is an in principle time-scale-dependent Haven raticP33SEU 
reflecting ion-ion correlations and 7 ~ 2 is a numerical factor reflecting the conductivity spectrum at onset of ac 
conduction, see Appendix A): 



Combining Eqs. Q and ([5| yields 



nq 2 (Ar 2 ft*)) , 



n A w 1 (Ar 2 ft*)) 
n 7 tf a 2 W 



If the mean-square displacement obeys (Ar 2 (t*)) ~ a 2 one has haw — ^mob(^ = V^*)> but unfortunately the 
quantity (Ar 2 (t*)) does not generally have this approximate value. In a simple models where all ions have similar 
jump rate, (Ar 2 (t*)) is indeed roughly a 2 times the fraction of ions that have jumped within time t*. It is not possible 
to model the universally observed strong frequency dispersion of the conductivity without assuming a wide spread of 
jump rates, however, and in such models like the random barrier model (RBM) considered below (Ar 2 (t*)) is much 
larger than a 2 . Generally, (Ar 2 (t*)) /a 2 gives an approximate upper limit for the fraction of ions that have moved in 
the time window t*. Ignoring the less significant factor jH, this implies that n mo b(^*) < %w- To summarize, only 
in models where all ions have similar jump rates does uaw gi ye an estimate of how many ions on average jump over 
the time interval of length t* . 

An alternative suggestion for obtaining information about the "number of mobile ions" is based on analyzing the 
electrode polarization regime of conductivity spectra for ion conductors placed between blocki ng electrod es 1=^^™ ^ 
However, theoretical analyses of the spectra are often based on Debye-Hiickel-type approaches l 42 l 43 l 44 l 45 l the applica- 
bility of which is far from obvious at high ion density. Thus while it is a potentially useful idea, more theoretical work 
is needed before observations of electrode effects may lead to safe conclusions regarding the number of mobile ions 
(see the next section that outlines the a simple approximate description) ; one still needs to specify the ti me scale that 
the number of mobile ions refers to. - Solid-state NMR methods such as mot ional narrowing experiment J 46 ! 47 ! 4 ^ and 
the analysis of multi-time correlation functions of the Larmor frequency] 49 * 50 * provide information about the number 
of ions moving on the time scale that these methods monitor (milliseconds to seconds). 

The question "what is the density of mobile ions?" is thus well defined only when it refers to a particular time scale. 
This is because according to standard ergodicity arguments, if the time scale is taken to infinity, all ions contribute 
equally and the density of mobile ions is the total ion density n. A natural choice of time scale is that characterizing 
the onset of ac conduction, the t* of the above equations. Choosing this time scale leads to a classification of ion 
conductors into two classes: Those for which n mo t>(£*) is comparable to the tota l ion density n: n mo b(t*) ~ n ("strong 
electrolyte case" 13 ), and those for which n mo b(t*) <C n ("weak electrolyte case"^^). The latter class includes solids 
where ion conduction proceeds by the vacancy mechanism (Sec. VII). 



III. WHAT CAN BE LEARNT FROM ELECTRODE POLARIZATION? 



As is well known, the ac conductivity cr(ui) is a complex function. Thus associated with the real part there is 
also an imaginary component; the latter determines the real part of the frequency-dependent permittivity. For the 
study of ion conduction in disordered solids the use of blocking or partially blocking metal electrodes is convenient. 
In this case, the high-frequency parts of ac conductivity and permittivity spectra are governed by ion movements 
in the bulk of the solid electrolyte, while the low- frequency part is governed by so-called "electrode polarization" 
effects, as shown in Fig. 3. Since the ions are blocked by the metal electrode, there is accumulation or depletion of 
ions near the electrodes, leading to the formation of space-charge layers. The voltage drops rapidly in these layers, 
which implies a huge electrical polarization of the material and a near-absence of electric field in the bulk sample at 
low frequencies. The build-up of electrical polarization and the drop of the electric field in the bulk are reflected in 
an increase of the ac permittivity and a decrease of the ac conductivity with decreasing frequency.^ For completely 
blocking electrodes <r(0) = 0, of course. - Whenever both ions and electrons conduct, a number of electrochemical 
techniques exist for ev aluating transference numbers of ions and electrons, including galvanic cells, polarization, and 
permeation techniques! 4 * 54 * 55 * 5 -^ 

Systematic experimental and theoretical studies of electrode polarization effects in electrolytes began in the 1950s 
in works was carried out by Macdonald 57 , Friauf 58 , Ilschner 59 , Beaumont 6 ^, and others. Their approaches were based 
on differential equations for the motion (diffusion and drift) of charge carriers under the influence of chemical and 
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FIG. 3: (a). At high frequencies the nearly constant loss (NCL) regime appears where the conductivity becomes almost 
proportional to frequency, sometimes even closer to a straightforward proportionality than shown here (data for a Lithium- 
Phosphate glass), (b) The electrode polarization effects on the real part of the conductivity and the real part of the dielectric 
constant at high temperature for a Na-Ca-phosphosilicate glass. 
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FIG. 4: Simplified electrical equivalent circuit describing the low- field ac conductivity and permittivity spectra of solid elec- 
trolytes between blocking (R c t = oo) or partially blocking electrodes. The right element describes the bulk properties, the left 
describes the space charge layer capacitance and charge transfer resistance (the frequency dispersion of the bulk conductivity 
may be taken into account by replacing the resistor Rb hy a, frequency-dependent impedance). 



electrical potential gradients. These equations were combined with the Poisson equation and linearized with respect 
to the electric field. Thereby, expressions for the ac conductivity and permittivity at low electric field strengths were 
derived. These are mean- field approaches in the sense that a mobile charge carrier interacts with the average field 
produced by the electrode and the other mobile carriers 61 . 

When charge carrier formation and recombination can be neglected and the sample thickness L is much larger than 
the space-charge layer thickness, the theoretical expressions can be approximately mapped onto the simple electrical 
equivalent circuit shown in Fig. 4 if the frequency dependence of the bulk conductivity is ignored. Ion transport in 
the bulk is described by the RbCb element. The space charge layers are described by a capacitance Cep , and, in 
case of a discharge of the mobile ions at the electrode, a (generally large) parallel charge transfer resistance R ct . The 
RciCep element acts in series with the RbCb element. In the cases Cep ^> Cb and R ct ^> Rb that usually apply, 
the equivalent circuit leads to the following expressions for the frequency-dependent conductance G{uj) (real part of 
the admittance) and capacitance C{uj) (imaginary part of the admittance): 



n(,,\ - (1/ Rct + Rb/ Filt) + u 2 R B C EP 

GH = M L = (I + Rb/^ + ^RbCep^ (7) 



and 
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= e ° € {UJ) L = CEP (l + R B /R^ + (.R B C EP y (8) 

with Cep/Cb = L/(2Ld) and A denoting the sample area, where the Debye length Ljj is defined by 

L D = — ~ 2 — ' W 

^mob^ 

From these expressions a number density of mobile ions n mo b can be calculated that is the density of mobile ions 
referring to the time scale for build up of electrode polarization, n mo b = ^mob(^ep) where r ep = RbCep- 

In the absence of ion discharge, i.e., when R ct — > oo , the equivalent circuit reduces to an RC element in series with 
a capacitor. The existence of a finite charge transfer resistance leads to the occurrence of a conductance plateau at 
low frequencies with plateau value G s given by 

Rb + R c t ^ ^ 

In addition, the static capacitance C s becomes slightly smaller than Cep- 

This mean-field approach should apply to materials with low n m ob? such as ionic defect crystals and diluted elec- 
trolyte solutions. Its applicability to disordered solids with high ion density is far from obvious. Nevertheless, quite a 
number of ac spectra of ion conducting glasses and polymers were traditionally analyzed and interpreted utilizing the 
above equations. Thereby, number densities of mobile ions were calculated and compared to the total ion content of 
the samples. For instance, Schiitt and Gerdes concluded that in alkali silicate and borosilicate glasses only between 
1 ppm and 100 ppm of the alkali ions are mobile^. Similar results were obtained by Tomozawa on silica glass with 
impurity ions^ and by Pitarch et al. from voltage-dependent measurements of on a sodium aluminosilicate glasJ^. 
Klein et al. carried out measurements on ionomers containing alkali ions and found ratios of mobile alkali ions to the 
total alkali ion content ranging from about 10 ppm to 500 ppm 44 . 

For a critical discussion of such experimental results and their interpretation, it is important to consider limitations 
of both experiment and theory. Regarding the experimental situation there are in particular two important points: (i) 
The roughness of the electrode/solid electrolyte interface is usually not taken into account. Especially in a frequency 
range where the length scale of the potential drop at the electrodes is comparable to the roughness of the interface, the 
roughness must have a considerable influence on the ac conductivity and permittivity, (ii) The surface-near regions 
of ion conductors often exhibit a chemical composition that is significantly different from the bulk. For instance, in 
ionic glasses, surface corrosion is initiated by an alkali-proton exchange. Such deviations from the bulk composition 
should have a strong influence on the ac spectra when the potential drop occurs very close to the surface, i.e., at high 
capacitance values close to the static capacitance plateau and in the static capacitance plateau regime. 

From a theoretical point of view, serious limitations of the applicability of mean-field approaches to disordered 
solids derive from: (i) The interactions between the ions are not taken into account; (ii) surface space charges in 
disordered solids may exist even without the application of an external electric field, due to ion exchange processes at 
the surface or due to an interaction of mobile ions with the metal electrode. Thus more sophisticated theories should 
take into account the possibility of an open-circuit potential difference between electrodes and solid electrolyte. 

In summary, considerable efforts in both experiment and theory is required in order to carry out measurements 
on well-defined electrode/electrolyte interfaces and to obtain a better theoretical understanding of what kind of 
information may be derived from electrode polarization effects. It is clearly worthwhile to pursue this direction of 
research, and it would also be worthwhile to look into what can be learned from electrode effects in the strong-field 
case where the electrode polarization becomes nonlinear. 



IV. WHAT CAUSES THE MIXED- ALKALI EFFECT? 



A prominent phenomenon occurring in ion-conducting glasses is the mixed-alkali effect (for reviews, se d 64 ' 65 ' 66 ' 67 ! ). 
This effect is the increase of the mobility activation energy of one type of ion when it is gradually replaced by a second 
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type of mobile ion. This leads to changes of the tracer diffusion coefficients Da,b(%) over several orders of magnitude 
at low temperatures, and to a minimum in the dc-conductivity a(x) ~ (1 — x)Da(x) +xDb(x), when the mixing ratio 
x of two mobile ions A and B is varied (for recent systematic experimental studies, see e.g. Refs. 68 69). 

Muc h progress has been made over th e last two decades for explaining the mixed-alkali 
effed 70 * 71 * 72 * 73 ^ Compelling evidence now exists that its origin is of structural 

character, associated with a mismatch effect!^ where sites in the glassy network favora ble for one type of mobile ion 
are unfavora ble for the other type of mobile ion. This evidence comes from EXAF^ESMl, NMRP and infrared 
spectroscop y 93 * 94 ^, x-ray and neutron scattering experiments in c ombin ation with reverse Monte-Carlo modeling^ 
and bond valence sum analysed, molecular orbital calculation d 96 * 97 1 , mol ecular dyna mics simulation d 74 * 98 * 99 ^ -! 
and theoretical work based on microscopic and semi-microscopic approache d 70 * 72 * 79 -^ 82 . In hopping systems, the 
mismatch effect can be modeled by site energies that are different for different types of mobile ions, i.e., a low-energy 
site for one type is a high-energy site for the other type. ____ 

Recently it was possible also to explain the peculiar behavior of the internal friction in mixed-alkali glassed^ 1 * 102 *. 
When a mixed-alkali glass fibre is twisted at a certain frequency, two mechanical loss peaks can be identified well 
below the glass-transition temperature: the single-alkali peak that with beginning replacement becomes smaller and 
moves to higher temperatures, and the mixed-alkali peak that at the same time becomes higher and moves to lower 
temperatures (for a review of experimental results, see^^). Based on general theoretical considerations it was shown 
that the mixed-alkali peak can be traced back to mutual exchanges of two types of mobile ions and the single-alkali 
peak to exchanges of the (majority) ion with vacancies. As a consequence, large mixed-alkali peaks are predicted 
for ion types with small mismatch where ion-ion exchange processes are more likely to occur. This agrees with 
experimental observations. Moreover, it could be shown that the occurrence of large mixed-alkali peaks at small 
mixing ratios can be understood if the fraction of empty sites is small. This gives independent evidence for the small 
fraction of empty sites found in theoretical argument s^ 4 ^ as well as in molecular dynamics simulation d 99 * 100 * 105 * 106 * 
(section VII). 

Despite this progress over the past years, there are still many issues awaiting experimental clarification and theoret- 
ical explanation. A point less addressed so far in the microscopic modelling is the behavior of the viscosity as reflected 
in a minimum of the glass-transition temperature upon mixing. This softening of the glass structure at intermediate 
mixing ratios may significantly influence ion transport properties. The mixed-alkali effect becomes weaker with total 
mobile ion content^ in agreement with theoretical expectations.^ However, a systematic theoretical study of this 
feature has not yet been undertaken. Overall, there is still no consistent theoretical account of all main signatures of 
the mixed- alkali effect. 

We finally note that a mixed-alkali effect also occurs in crys tals with s tructure of (3- and /^"-alumina type, where the 
ion motion is confined to two-dimensional conduction plane s! 107 * 108 * 109 * A quantitative theory has been developed for 
this based on the wealth of structural information available! 110 * 111 * The key point is that A and B ions have different 
preference to become part of mobile defects, and this preference is caused by a different interaction of the ions with 
the local environment. In this respect the origin of the mixed-alkali effect in crystals has similarities to that in glasses. 
However, different from the host network in glasses, the host lattice in the crystalline systems is almost unaffected by 
the mixing of the two types of mobile ions. 

V. WHAT IS THE ORIGIN OF TIME-TEMPERATURE SUPERPOSITION? 

Different suggestions were made in the past to characterize the similar ac responses observed for different types of 
ion conductors in frequency regimes not exceeding ~100 MHz. The simplest description is the power-law frequency 
dependence proposed by Jonscher (Eq. ([3]) ) l 112 l 1= E 51 The power-law description is not accurate, however, because the 
exponent must generally increase somewhat with frequency in order to fit experiment properly, and also because 
the asymp totic low-frequency behavior is inconsistent with experiment that imply cr'(u) — cr(0) oc uo for co — > 
(cj = 2ttv)W^^ A more general appro ach is to c onsider the scaling associated with time-temperature superposition 
(TTS) for any particular ion conducto J 115 * 116 * ^ ^. The scaling ansatz reads 

a(u;,T) = a dc (T)f[ou/ou*(T)] . (12) 

Here f(u) is the so-called scaling function and ou* the previously defined angular frequency marking onset of ac 
conduction. Any solid that obeys TTS is, equivalently, referred to as obeying scaling. As an example Fig. 5(b) 
illustrate s how the spec tra of Fig. 2(b) scale to a common so-called master curve. 

BartonpU NakijamapHl and NamikawsP^ long ago verified that for many ion- (and some electron-) conducting 
disordered solids 
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FIG. 5: Time-temperature superposition (TTS). (a) Sketch of how the three spectra of Fig. 2(a) collapse to a single master 
curve when suitably scaled. Whenever this is possible, the solid is said to obey TTS. Most disordered ion conductors, including 
all single-ion conducting glasses, obey TTS. (b) TTS demonstrated for the Lithium-Phosphate glass ac data of Fig. 2(b). 
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where p is a constant of order of unity and Ae is the dielectric strength, i.e., the difference b etw een static and high- 



frequency dimensionless dielectric constants. Equation (13) is known as the BNN relation By considering the 



low-frequenc y exp ansion of the conductivity a connection can be made between the scaling behavior Eq. (12) and the 
BNN relation!^] The argument assumes analyticity of the scaling function f(u) for small u, which is in fact necessary 
in order to have a well-defined dielectric strengthP^ One has cr(uj)/cr& c ~ 1 + iKu/u* for u — > with the constant 
K being real. Accordingly, one obtains e(uo) — = <j{uj) / {ie^uS) ~ crdc/faow) + Ka& c j e§uj* for uo — > which implies 
Ae = K(j& c / e^ou* . Thus TTS and analyticity imply the BNN relation - but do not mathematically garantee that 

To the best of our knowledge, TTS applies for all single-ion conducting glasses and crystals with structural disorder. 
This remarkable fact suggests that disorder is intimately linked to TTS. In crys tals with s tructural disorder, such as 
RbAg4l5 and /3-alumina, different types of ion sites exist with different energies! 124 * 125 * 12 ^ In addition, the interionic 
Coulomb interactions cause a significant spread in the potential energies of the ions. In glasses, the disorder of the 
glass matrix l eads to a broad distribution of ion site energies and barrier heights and thus to a broad distribution 
of jump 

ratesPZEB] This may explain why, even in single-modified glasses (i.e., with only one type of ion) with 

low numb er ion density and corresponding weak interionic Coulomb interactions, violations of TTS have not been 
observed.™^ 

In contrast, crystals with low concentrations of point defects routinely show TT S violations. Examples are materials 
with intrinsic Frenkel or Schottky defects, such as alkali and silver halidesP^ In these materials, the interactions 
between the small number of defects are weak and the defects are partly bound to counter charges. Therefore, on 
short time scales the defects carry out localized movements close to the counter charges. These localized movements 
are not correlated to the long-range ion transport, and consequently the conductivity spectra do not obey TTS. 

Violations of TTS are also found in materials with more than o ne type of mobile ion. Examples are mixed- 
alkali glasses^], as for instance 2 Ca(N0 3 ) 2 • 3 KN0 3 (CKN) melts liMSU Below the glass transition temperature 
(T g = 333K) CKN is believed to be a pure K + ion conductor and it obeys TTS, but at higher temperatures Ca 2+ 
ions most likely also contribute significantly to the conductivity (above 375K CKN again obeys TTSM^). Other 
examples are some polymer electrolytes above their glass-transition temperature where the polymer chains carry out 
segmental movements. Here different types of movements with different characteristic length scales contribute to the 
conductivity spectra, which generally results in TTS deviations p2U 



The simplest model exhibiting the scaling properties Eqs. (12) and (13) is perhaps the random barrier model 
(RBM), see, e.g., Ref. 36 for a review. In this model hopping of a single particle on a lattice with identical site 
energies is considered, where the energy barriers for jumps between neighboring sites are randomly drawn from a 
smooth probability distribution. The particles must overcome a critical "percolation" energy barrier E c to exhibit 
long-range motion (a review of percolation theory with particular emphasis on ion diffusion was given by Bunde 
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FIG. 6: Approximate ac universality and deviations from it. Figure 6(a) shows the RBM prediction for the scaled real part of 
the ac conductivity (full curve) compared to data for two typical ion conducting glasses, the 24°C data of Fig. 2(b) and the 
Sodium-Borate data of Fig. 6(b). The fit is good, but not perfect. Figure 6(b) shows that ac universality does not include 
mixed-alkali glasses (blue); the red curve for a Sodium-Borate glass represents approximate ac universality, compare Fig. 6(a). 



and HavlirPl). The time to overcome the percolation barrier, t c cx exp(E c /kBT), determines the characteristic 
frequency marking onset of ac conduction: cu* ~ t" 1 .^ The percolation energy barrier - acting as a bottleneck - 
also determines the dc conductivity temperature dependence. Thus percolation explains why a wide distribution of 
barriers nevertheless results in an Arrhenius dc conductivity (which is observed for most disordered ion-conducting 
solids). Incidentally, the BNN relation's rough proportionality, <r(0) ex cj*, also follows from percolation determining 
the conduction properties. 

The scaling function of the RBM is universal in the "extreme disorder limit" where the jump rates vary over several 
decades; this limit is approached as temperature is lowered. Universality means that the ac response in scaled units 
becomes independent of both temperature and activation energy probability distribution. This was shown by extensive 
computer simulations invo lving barriers distributed a ccor ding to a Gaussian, an exponential, an inverse power law, 
a Cauchy distribution, etc! 36 * 137 J It was recently showrP^that if a = a(cj)/ad c and a; is a suitably scaled frequency, 
except at low frequencies where the conductivity approaches the dc level, the universal RBM ac conductivity is to a 
good approximation given by the equation 



/ \ 2/3 

. (14) 

This expression implies an approximate power-law frequency dependence of the ac conductivity with an exponent 
that slowly converges to unity at very high frequencies - not simply an exponent of 2/3 as one might naively guess. 
A very accurate representation of the RBM universal ac conductivity is given in Appendix B . 

The RBM scaling function is usually close to, but rarely identical to those of experiments! 36 * 79 ^ 3 -^ As an example 
Fig. 6(a) shows the RBM universal ac conductivity (full curve) with the Lithium-Phosphate data of Fig. 2(b) at the 
lowest temperature (24°C) and the Sodium-Borate data of Fig. 6(b), where both data sets were empirically scaled 
on the frequency axis. Figure 6(b) shows data also for a mixed-alkali glass; these data deviate significantly from 
the approximate ac universality represented by the red curve. It appears that, on the one hand, the RBM captures 
the essential features of the ion dynamics for single-ion conducting disordered solids, and that, on the other hand, 
deviations from the RBM universal ac conductivity may provide important information about specific features of the 
solid in question 38 . Thus the RBM may be regarded as the "ideal gas" model for ac conduction in disordered solids. 

In the RBM the dispersive tra nsport properties are governed by strong disorder, forcing the ions to explore percola- 
tion paths for long-range motionP 40 * 141 ! Macroscopic alternatives to the R BM (but with similar basic physics 36 ), which 
apply if the sample has microstructure, have also been studied P 42 * 143 * In most cases except that of nanocrystalline 
materials, however, the disorder is believed to be on the atomic scale. 

Single-particle models like the RBM or its generalizations are simple and attractive for understanding the origin of 
TTS, but there are a number of challenges to this approach as well as open questions that one must keep in mind: 
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(%) For glasses the stoichiometry can be varied to a large extent. Related to changes in composition there are 
changes in activation energies, as for example a lowering of the activation energy with increasing mobile ion 
content or an increase of activation energy when one type of mobile ion is successively replaced by another type 
(the mixed-alkali effect, Sec. IV). These effects are not dealt with in the RBM unless the model is modified in an 
ad hoc manner to allow for significant changes in the barrier distribution and mismatch effects with respect to 
different ion types. The simplest models accounting for these effects are hopping systems with site exclusion) 79 * 80 ! 
i.e., where there can be at most one io n at each site. Interestingly, such "Fermionic" hopping systems with site 
energy disorder often obey TTS! 144 * 145 ! Moreover, calculations for the corresponding single-particle systems yield 
scaling functions^!! that are close to the RBM universal scaling function. A conclusive picture of the scaling 
properties of these type of models remains to be established, however. 

(ii) Recent mol ecula r dynamics simulation d 99 * 100 * 105 * 106 * 14 ^ and theories for the internal friction behavior in mixed- 
alkali glasses^J show that often only few of the potential ion sites are vacant (typically 1 — 3%). This is expected 
on gen eral grounds, since during the cooling process a glass tends to a state of low free energy, thus with few 
defect sP^ It would be interesting to investigate whether hopping models with a low concentration of vacant 
sites generally obey TTS. 

(Hi) The Coulomb interactions between ions can be estimated from their mean distance R oc n _1//3 , where n is the 
number density of ions. At room temperature, typical plasma parameters e 2 [( 4:'Ke 00 Rk&T) are in the range 
30-80. In view of the confined geometry of the diffusion (percolation) paths 148 the local interactions may be 
even stronger. Hence it is important to clarify whether hopping models with Coulomb interaction obey TTS 
and, if so, how the scaling function is affected by the interactions (for a general overview of Coulomb interactions 
effects on dispersive transport properties, see e.g. Ref. 1123)) . Early studies of Coulomb interaction effects in 
hopping models with percolative disordep^ showed that Coulomb interactions give rise to a strong conductivity 
dispersion, but TTS was not observed. This might be due to the fact that in these early simulations temperature 
was not low enough. Another reason could be that, as in the RBM, a smooth and broad distribution of barrier 
or site energies is required for scaling. Indeed, simulation studies of many-particle hopping in the RBM with 
Coulomb interactions show agreement with the scaling behavior for low and moderate particle concentrations 
in the limit of low temperatures Overall, however, the problem is far from being settled; in particular if 
one takes into account that the fraction of empty sites should be small and that critical tests for other types 
of structural disorder such as site energy disorder, have not yet been performed. Due to the long-range nature 
of the Coulomb force, one could argue that its contribution to the energy landscape (sites and saddle points) 
provides an overall mean-field contribution. This hypothesis should be tested by further simulations. 

(iv) Most studies of the RBM and other hopping models focused on site an d/or ba rrier energies varying randomly 
without spatial correlation (a notable exception is the counterion modet^ 4 ^ 5 -^) . If there are significant spatial 
correlations - thus introducing a further length scale into the problem - this may well lead to a breakdown of 
TTS. 



VI. WHAT CAUSES THE NEARLY CONSTANT LOSS? 

At high frequencies and/or low temperatures conductivity spectra approach a regime with nearly linear frequency 
dependence when plotted in the usual double- log plot: (r'(uS) oc Auf 1 (n = 1). The proportionality constant A is only 
weakly temperature dependent. This is referred to as the "nearly constant loss" (NCL) regime since it corresponds to 
an almost frequency-independent dielectric loss (Fig. 3). This behavior is ubiquitously observed in a wide variety of 
solids including glassy, crystalline, and molten ion conductors, independent of specific chemical and physical structures 
- for an overview on experimental results, see e.g. Refs. I151|1521 

There are different possible origins of the NCL. One possibility is that NCL reflects the still not fully understood 
low-energy excitations present in all disordered materials. In the quantum-mechanical tunnelling model these excita- 
tions account for the anomalous low-temperature features of heat capacity and sound-wave absorption!^ 5 - 3 -* At higher 
temperatures the low-energy excitations give rise to relaxations of the system over an ene rgy barrier separating two 
different energy minima, described by the asymmetric double- well potential (ADWP) modelP^On a micros copic level 
this could correspond to cooperative "jellyfish-type" movements of groups of atomic species in the material! 135 1 155 1 156 ^ 
If correct, this type of dynamic process should be a feature of all disordered materials, including materials without 
ions. 

A more recent interpretation suggests that localized hopping movements of ions within fairly small clusters of sites 
contribute to the NCL in disordered ion conductors! 157 * 15 ^ In this interpretation the NCL is merely the extension 
to higher frequencies of the dispersive conductivity. In fact, any hopping model with sufficient disorder gives rise to 
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Vacancy-like transport mechanism 




FIG. 7: Three vacancies jump to the left, each proceeding via a sequence of ion jumps to the right. First one ion jumps into 
the vacancy, then the next ion jumps into the new vacancy, etc. 

such a regime, since on short time scales hopping models always correspond to ADWP-type models. In the RBM, 
for instance, ion jumps over limited ranges lead to an NCL regime at high frequency, a region that extends to lower 
frequencies as temperature is lowered. There are two experimental observations favoring the second interpretation: 
(i) The magnitude of the NCL increases with increasing ion concentration^^ (this also applies in the ADWP model 
if the defect centers are somehow associated with the mobile ions); (ii) at large ion concentration and temperatures 
above 100 K, the scaling properties of the NCL contribution to the conductivity spectra are identical to t he sca ling 
properties found at lower frequencies where the dispersive conductivity passes over to the dc conductivity ! 1 57 * 1 ^ On 
the other hand, experiments carried out by one of the present authors suggest that the low-temperature (T < 80 
K) NCL in glasses with few ions is due to ADWP-type relaxation of the glass network This i ndic ates that both 
hopping movements of the ions and ADWP-type relaxations in the material contribute to the NCLP^Sl Which of these 
dynamic processes dominates depends on composition and temperature. 

Again, one may wonder whether it is permissible to neglect interactions, which can be modelled by dipolar forces 
as regards the short-time dynamics with only local movements of the mobile ions close to some counterions. Monte 
Carlo studies and analytical calculations have shown that the energetic disorder associated with the local electric field 
distributions of spatially randomly distributed dipoles gives rise to an NCL contribution at very low temperatures 
within an effe ctive one-particle descri ption, whereas at higher temperatures such behavior can occur due to many- 
particle effects! 149 ' 150 ' 16 ^ 162 ' 163 ' 164 ' 165 ' 

Finally, it has been suggested that the NCL is caused by vibrational movements of the mobile ions in strong ly anhar- 
monic potentials or from ion hopping in a slowly varying cage potential defined by neighboring mobile ions! 166 l 167 l 168 l 
These views focus on the very high-frequency NCL. Indeed, at frequencies in the THz range the ac conductivity joins 
into the v ibra tional absorption seen in far infrared spectroscopy associated with the quasi-vibrational motion of the 
mobile ion 151 Unfortunately, the connection between the NCL and the vibrational modes is poorly investigated: A 
data gap from the GHz to the THz regions exists because measurements of the dielectric loss are here particularly 
challenging. More focused studies in this frequency window are needed to elucidate the connection between vibra- 
tional and librational (anharmonic) motion, as well as to better characterize the precise frequency dependence of 
the NCL cond uctivity, i.e., is it exactly linear (n = 1), slightly sub-linear (n < l,n = 1) or slightly super-linear 

(n>i,n^i) < PMnnn2] 

VII. WHAT IS THE ION TRANSPORT MECHANISM? 

We now turn briefl y to the most fundamental question relat ing to ion conduction in disordered solids: What is the 
transport mechanism^ 1 * 11 * 12 * 14 * 17 * 19 * 173 * 174 * 175 * 176 * 177 * 178 * 179 * 180 -! As mentioned already, there is evidence that in many 
cases only few of the sites available for ions are actually vacant. Figure 7 illustrates the vacancy mechanism. Ion 
conduction in disordered solids does not proceed via the well-understood vacancy mechanism of ionic crystals with 
few vacancies In crystals the vacancy concentration is strongly temperature dependent due to thermally activated 
defect formation; in glasses the concentration of empty sites is determined by the history of glass formation via the 
cooling rate, and the number of vacancies is frozen in at the glass transitions^. Moreover, the vacancy concentration 
in glasses is believed to be significantly larger than in crystals. This makes the conduction mechanism in glasses 
much more complex, also because vacancy-vacancy interactions generally cannot be ignored. If such interactions are 
nevertheless not important, one may regard a vacancy as a charge carrier, e.g., in the RBM. In this case, the results 
of the RBM pertaining to the scaling features of conductivity spectra apply unaltered. 

Evidence for the significance of many-particle effects als o comes from molecular dynamics simulations and from 
measurements of the Haven ratio. Simulation d 105 * 182 * 183 * 184 * suggest that hopping occurs in a cooperative way, where 
one ion jump triggers jumps of other ions or hopping events occur collectively by involving several ions. A directional 
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correlation of jumps of different ions is indicated by the measured Haven ratios that are generally smaller than 
unityf^ corresponding to positive cross-correlations in the current autocorrelation function. If heterogeneities in the 
host networ k confi ne the ion motion to channel-like structures, such correlations can be expected to be even more 
pronounced! 148 * 185 ^ Haven ratios smaller than unity have also been found in Monte Carlo simulations of models with 
Coulomb interactions 

There is a need for more systematic investigations of the role of many-particle effects. To uncover the ion transport 
mechanism (s) model predictions should be compared to other experimental observables than conductivity, such as 
NMR, spin-lattice rel axation, mechanical relaxation, tracer diffusion, quasi-elastic neutron scattering and multiple 
spin-echo experiment J 186 l 187 l 188 [ These methods probe different correlation functions, so checking model predictions 
against them obviously provides more severe tests than just, e.g., comparing predicted ac conductivity spectra to 
experiment. Such lines of inquiry, however, have so far only been undertaken in few instances; for example spin-lattice 
relaxation spectra were investigated in Refs. 35 50 189 , and multiple spin-echo experiments in Ref. 11901 - The 
ion conduction mechanism clearly deserves a review on its own summarizing the latest developments. We have here 
mainly announced the problem and left out detailed considerations of, e.g., the role played by dynamic heterogeneities 
and possible dynamic channels for conduction pathway d 184 * 191 *. Hopefully answering the other questions of the present 
review will provide valuable input into revealing the ion conduction mechanism. 



VIII. WHAT IS THE ROLE OF DIMENSIONALITY? 



The subdiffusive ion dynamics on transient time scales found ubiquitously in disordered materials is sometimes 
attributed to a fractal geometry of the conduction paths. A classic example of subdiffusion is particle dynamics 
occurring on a percolation cluster^ 6 * 192 *. In this case the cluster is fractal on length scales below a correlation length 
£, which diverges when approaching the percolation threshold. At criticality, dangling ends and loops occur on all 
lengths scales causing the mean-square displacement to increase as a power law with an exponent smaller than one. 
Close to the percolation threshold, £ is finite and subdiffusive behavior is observed only in an intermediate time 
regime, where the mean-square displacement is larger than microscopic length scales and smaller than £ 2 . For times 
where the mean-square displacement exceeds £ 2 , the diffusion eventually becomes normal. This examples suggests 
that the effective dimensionality of the conduction pathways may play an important role for the subdiffusive behavior, 
although for conduction pathways containing loops the fractal dimension of the pathway structure and the embedding 
Euclidian dimension of the material are generally not sufficient to determine the power law exponent in the subdiffusive 
time regime (see the discussion in chapter 3 of ref. 136). Hence the question arises to what extent, if any, does the 
dimensionality of the conduction space influence the subdiffusive motion? 

A possible scenario for qualitatively understanding the origin of the dimensionality dependence is the following. 
With decreasing dimensionality, the average distance between the highest barriers (percolation barriers) on the con- 
duction pathways becomes larger. Between these percolation barriers, the ions perform back- and- forth motion. An 
increasing spatial extent of this back-and-forth motion leads to a larger dielectric relaxation strength, implying a more 
gradual transition from dc conductivity to dispersive conductivity in the conductivity spectrum. 

Few studies of this question exist, but there is some evidence that dimensionality does influence the shape of the 
ac conductivity master curves (the scaling functions of Eq. ([12]) ) . This can be seen in 2D crystals, like sodium 
(3 — AI2O3, and in ID crystals, like hollandite: The transition from dc conductivity to dispersive conductivity becomes 
more gradual with decreasing dimension!^' This variation in the shape, as characterized by an effective exponent of 
Eq. (14), is shown in Fig. 8. This sensitivity to dimensionality is also evident in the RBM for which the shape of 
the conductivity spectra are similarly altered by changing the dimensiona lity, i.e ., in two dimensions the conductivity 
increases somewhat less steeply with frequency than in three dimensions ! 194 * 19 ^ 

The conducting pathways of the two above-mentioned crystals (sodium f3— AI2O3 and hollandite) have a well-defined 
dimensionality. Is there any other evidence for modifications of the correlated motion which might be connected to 
the dimensionality (localized) of the conduction space in an amorphous solid? While mos t disordered materials as 
mentioned show very similar shapes of the ac conductivity curves, some researcher jj 196 l 197 Il£8] have observed systematic 
changes that appear to arise from changes in the local environment of the ions. In studies of alkali- germanate glasses 
and alkali-borate glasses of varying ion content, for instance, subtle changes i n the shape of the conductivity spectra 
were seen to correlate to known anomalies in the glass-transition temperature! 37 * 38 -* The glass-transition temperature 
passes through a maximum as a result of how added modifier ions initially polymerize, but later depolymerize, 
the oxide network. Consequently, the average oxygen coordination of the ion's charge-compensating sites changes 
with ion concentration, resulting in modifications of the local ion environment which could mimic changes to the 
dimensionality of the conductio n spa ce In a similar way, the mixing of two ion specie s (s ay Li and Na) modifies 
the local environment of the ionp^* and studies of ac conductivity of mixed-alkali glasses 196 show a distinct change 
compared to that found in either single alkali end- member composition alone (compare Fig. 6(b)). Finally, in studies 
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FIG. 8: Dimensionality dependence of the conductivity exponent n (Eq. (pi)) for glasses, planar crystals and channel crystals. 



of metaphosphate glasses (whose oxide structure is highly polymeric) systematic changes in the correlated motion 
occurred in conjunction with variations in the cation size relative to the free volume*™. There it was posited that an 
effective local dimension of conduction space might rule the correlated motion. - More studies along these lines are 
clearly warranted, but it appears that the effective dimensionality is an important parameter. 



IX. CONCLUDING REMARKS 



Science progresses by asking questions. It is our hope that this paper will stimulate to works focusing on basic 
understanding, eventually leading to a physical picture and quantitative model(s) of ion conduction in disordered 
solids that are as good as those of ionic and electronic conduction in crystals. It is a reasonable working hypothesis 
that ion conduction in disordered solids can be described in terms of a fairly simple generally applicable model, but 
only the future can tell whether this hope is realized. 
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APPENDIX A: RELATION BETWEEN THE LONG-TIME MEAN-SQUARE DISPLACEMENT AND 
THE LOW-FREQUENCY BEHAVIOR OF THE AC CONDUCTIVITY 

The Kubo formula in dimension d (where s = ico + e is the "Laplace frequency" and e > is eventually taken to 
zero) reads 

^) = ^i *d(0) ■!(*))«-. (Al) 
Here V is the sample volume and the total current I(t) is defined by summing over all N ions: 



N 

I(t) = qJ2vj(t). 



(A2) 
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Introducing the collective displacement over a time interval of length t, 



we have 



N 1 r t 

AR(t) = ^ Ar,-(t) = - y dr'I(r') . 



(A3) 



q 2 (AR 2 (t)) 

2 d(AB?(t)) 
1 dt 



= f dt! f dt 2 (I(ii) • I(i 2 )> =2 f dr(t- r)(I(0) • I(r)> , 
Jo Jo Jo 

= 2jfdT<I(0).I(r)>, / 2(A ^ (t)) = 2(1(0) ■!(*)). 



Accordingly, after a partial integration where t he b oundary term disappears because AR 2 (t) ~ t 2 for £ — > (reflecting 
the short-time socalled ballistic motion), Eq. (Al) takes the form 



a(s) 



2dVk B T 



f 

Jo 



eft 



C 



Jo 



dtf(t), 



(A4) 



where C = nq 2 /k B T and /(t) = (AR 2 (t))/2diV. 
We now make the ansatz 



(») ~ ^dc + -4s" 



(A5) 



with A real and < a < 1. This is the well-known Jonscher ansatz^- 2 * 113 * ana lytically continued to complex frequencies, 
because for the real part of the frequency-dependent conductivity Eq. (A5) implies ct'(uj) ~ 0"dc[l + (ui/uj*) a ] with 



(Jdc(u*)- a = Acos(mr/2). 



(A6) 



The analyticity requirement implies that at very low frequencies one must eventually have a = 1, but as an effective 
description of the regime of onset of ac conduction this ansatz may still be used. It follows that 



cr(s) - a dc 



/ dt Cf(t) 
Jo 1 



As c 



0. 



Thus via a Tauberian theorem one concludes that 



(A7) 



Cf(t) - <7 dc 



t — > oo . 



r(i-a) ' 

Since asymptotic expansions can be integrated term by term and /(0) = 0, we obtain 



(A8) 



Cf(t)~a dc t 



A 



T(2-a) 



l-a 



t — > OO 



or, if D(j = <Jdc/C — (^bT '/nq 2 )cr ( ±c — lim^ 00 (AR 2 (t))/6A/'£ is a collective diffusion coefficient, 



(A9) 



(AR 2 (rQ) 
2dN 



(w* t) 1 



- cos(a7r/2)r(2 - a) 



t — > oo . 



(A10) 



If we introduce the time scale t* = 1/z/* = 2tt/uj* corresponding to the above-defined crossover frequency uj* where 
a'(uj*) = 2<Jdc, and if we assume that the Jonscher ansatz is a good description near the crossover so that the 
asymptotic expression applies, Eq. ( |A1Q| ) implies 
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V (AR 2 (t*)} 

c = 2dk B T 7 N V (AH) 



with 



7 = 1 + 7 — x ; / — \ — 7 ^ • (A12) 

; (27r) a cos(a7r/2)r(2-a) v } 

The factor 7 is roughly two for a < 0.82, but diverges for a — > 1. If we further replace the collective mean-square 
displacement (AR 2 (t*)) by the single-particle mean-square displacement (Ar 2 (£*)), using the definition of the time- 
dependent Haven ratio 

1 l E j¥fc (Ar,(r)-Ar fc (r)) 
#(f) TV (Ar 2 (i*)) 1 ; 

that can b e appro ximated by the Haven ratio H in the dc-limit that is accessible via radioactive tracer 

experiments,™EIl i.e., (AR 2 (i*))/iV = (Ar 2 (t*))/H(t*) ~ (Ar 2 (t*)}/H, we arrive at Eq. of the main text 
(where d = 3): 

V (Ar 2 (t*)) 

adc= 2& B T 7 tf 17 * (AM) 

APPENDIX B: ACCURATE REPRESENTATION OF THE RANDOM BARRIER MODEL UNIVERSAL 

AC CONDUCTIVITY 



The expression Eq. ( 14 ) gives a good overall fit to the universal ac conductivity of the random barrier model (RBM) 
arising in the extreme disorder limit, but in the range of frequencies where the conductivity approaches the dc level 
(uo < uo*) there are significant deviations for both the real and imaginary parts of the conductivity. Thus for the 
imaginary part cf"{<jS) Eq. |l4| predicts that a" oc a; 2 / 3 for Co — > 0. This contradicts experiment, analyticity as well as 
RBM computer simulations"™ that all imply a" oc Co for uj — » 0. An accurate analytical representation of the RBM 
universal ac conductivity with the correct low-frequency behavior is giverP^ by the equation 

Ina= y ( k l + - y j . (Bl) 
In this expression frequency is scaled such that a = 1 + i& for Co 0, i.e., the frequency scaling is different from that 



of Eq. (14). This equation is easily solved numerically for both real and imaginary parts as functions of frequency. 



Alternatively, numerical solutions to it - as well as to Eq. (14) - are available as ASCII files (see reference 22 of Ref. 
fT38l) . 
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